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foreword 

The material presented in this final report was performed under &-ant 
No,. NAG-1-301 entitled "Component Mode Synthesis and Large Deflection Vi- 
brations of Complex Structures." This report sutraarizas the research re- 
guilts on modal synthesis and nonlinear forced vibrations of beams. The 
study was performed at the NASA/Langley Research Center during the period 
from Novenber 1, 1982 to October 31, 1983. The work was monitored under the 
supervision of Joseph E. Walz and Dr. Jerrold M. Housner, Structural 

Oynamcs Branch, Structures and Dynsnics Division, NASA/Langley Research 
Center. 
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COMPONENT MODE SYNTHESIS AND LARGE DEFLECTION VIBRATIONS 
OF Ca^PLEX STRUCTURES 


By 

Chuh Itel* 

Part 1. Dynamic Analysis of Large Complex Structures 
Using Component tode ffethods in NASTRAN 

The canplexity of aerospace structures has been increased enormously 
during the past decade. A new challenge has confronted the structural 
dynamists by the proposed space station to be in service by the year 1990. 

It will be an evolving structure (ref. 1), and it will not be possible for 
it to be ground tested because the final configuration may not be known when 
the first component is put into space. The component mode method, there- 
fore, may be onployed for the dynanic analyses for deteraining frequency, 
mode shape and transient response of such a large structure system in space. 

The NASTRAN computer program, a structural analysis tool widely used in 
the aerospace industry, contains a modal synthesis capability. Other than 
the nine-bay truss structural problen presented in the NASTRAN demonstration 
manual, little is publicly known about its capabilities. Preliminary 
assessment of the accuracy of the NASTRAN modal synthesis analysis is accan- 
plished by making a comparison of the NASTRAN modal synthesis with full 
structure NASTRAN and nine other modal synthesis results (ref. 2) using the 
nine-bay truss shown in Figure 1. Figures 2-4 show the relative accuracy 
obtained using the various modal synthesis procedures. The limited study 
indicates that the fixed-interface method in NASTRAN, fixed-interface 

‘^l^fon'afe TFoTessoF, Dspartaent of tiachanical Engineering and Mechanics, 

Old Dominion University, Norfolk Virginia 23508. 
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Figure 1. Truss Model 
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Comparison of Methods with Frequency Error of 0.1% 











Figure 3. Comparison of Methods with Frequency Error of 0.5%. 
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method of Kurty (ref. 3) and the equivalent Craig-Bampton method (ref. 4) 
produce the most accurate results for a given nunber of degrees-of-freedon. 
Slightly better accuracy was achiTsved by the procedures introduced by 
Benfield and Hruda (ref. 5). But these latter methods (ref. 5) suffer the 
disadvantage that the modes of one substructure are not independent of the 
modes of other substructures, f-fere detailed results are documented in the 
progress report entitled, "NASTRAN Modal Synthesis Capability (ref. 6)." 

A NASTRAN coinponent mode transient response analysis was also performed 
on the free-free truss structure. A concentrated force was applied at grid 
point 42 of component B for 0.12 seconds and then removed. Tables 1 and 2 
give the displacements and stresses and are ccmspared with the full structure 
NASTRAN results. It danonstrates that excellent transient response can be 
obtained using component modal synthesis, (tetail results, KiAP alters, sol- 
ution sequences and conputer CPU time can be found in reference 6. 
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Table I. Transient Response and Percent Error in Displaceaent 


O o 
-n ^ 

o a 

O :■;» 


Full Truss i S ^ubstr. { S I Full Trees 


g. Se&str. 


F*EI Full Truss 
F 


B Fubstr. ' F-BI 
? 


O.OOiS 

0.0030 

0.0045 

o.oo&a 

0.0073 
0 .0090 
G.0S05 

a .0120 


2 @ 

0.0 

0. 473^243 
2.Q709Q6 
4.794056 
8.662573 
13.65921 
19.79589 
27.07146 
35.47588 


41 

0.0 

0.4785015 
2.070847 
4.7938S2 
S. 662563 
13.65901 
I9.795S7 
27.07133 
35.47573 


0.0 

0.0 

0,0 

0.0 

0.0 

0.0 . 

0.0 

0.02566 

0.4S9332 

©.43S3054 

O.OS533 

6. 4736243 

0.476S01S 

0.02366 

O.C02S3 • 

2.oassi5 

2s0i9§3i 

-s.esiiQ 

I.STOSSS 

2.87SS4? 

O.0023S 

0,00363 

4.ai371i . 

4.613SS3 

0.00334 

4, 79.4056 

4.7§3S32 

0.03363 

0.00012 

8.6S2S71 

6. 663060 

-0.00218 

@.iS23?3 

S. 662563 

0.SC3I2 

0.00146 

13.67006 

I3.677S1 

0.00133 

13.63921 

13.63901 

0.00146 

O.OOSIO 

19.61603 

19.SI6I9 

-6.C0S5Q 

19.7SS89 

I9.79537 

0.03310 

0.Q3048 

t7.Q9IIS 

27,09121 

-o.oosit 

27.07146 

27.07133 

0.0S34S 

0.00042 

35.43301 

3S. 43476 

0.00070 

35.475S3 

35.47573 

O*GS'042 











PART 2, Large [teflectiosi Vibrations of Ba.sns using Finite Elesaent ^tethods 

Since space structures ^ill be mostly large, lightweight and flexible, 
large deflection analysis methods are urgently needed to study the v1” 
bratory responses of complex structures. 

A finite element method has been developed for nonlinear vibrations of 
bean structures subjected to harmonic excitation. Longitudinal deformation 
and inertia are both included in the formulation. A harmonic force matrix 
[h] was developed for a beam elonent for nonlinear oscillations under uni- 
form harmonic excitation. Formulation of the harmonic force matrix is based 
on the mathematical basis (ref. 7) that the simple harmonic force Cosut 
is simply the first order approximate solution of the simple elliptic fore" 
ing function BA cn(px,k). Also the well known perturbation solution 



] of a Duffing system + q ■}• 3q^ = Cc»swt is the first order 

approximate solution of the simple elliptic response q = A cn(px,k). Dari- 

A 

i 

I vation of the element harmonic force and nonlinear stiffness matrices are 

given in detail in progress report entitled, “Finite Element Analysis of 
Nonlinear Free and Forced Vibrations of Beisns (ref. 8)." 

Table 3 shows the finite elenent free vibration results with and with- 

I out considering effects of longitudinal deformation and inertia (ELDI). It 
clearly demonstrates the remarkable agreement betwen the present finite 
elonent with ELDI and Rayleigh-’Ritz solutions. 

Table 4 shows the frequency ratios for a simply supported and a clamped 
bean (L/R = 100) subjected to an uniform harmonic force of = 2.0 = 

1322 Ib/in.) and 1.0 (3277 lb/1n.), respectively. It demonstrates the 
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Function 
Solutloa 
(ref. 9) 

First 

Iteration 

Final 

Result 

lies 
Solution 
(ref. 10) 

First Pinal 

Iteration Result 



1.0 

1.0S92 

l.OSfS 

1.0888 

1.0607 

1.0613 

1 0613(3)^ 

2.0 

1.3178 

1.3203 

1.3119 

1.2246 

1.2270 

1.2269(4) 

3.0 

I.S2S7 

1.6295 

1.6022 

1.4573 

1.4620 

1.4617(4) 

4.0 

1.5760 

1.9761 

1.9216 

1.73Q9 

1.7383 

1.7375(6) 

.5.0 

2.3501 

2.33% 

2.2544 

2.0209 

2.0393 

2.0370(7) 


a. Effect© of loQgitudiiaal dsforaatioa aM iaoreia. 

b. I^i£sib@r in bracket© denote© the nnaber of iteration© to gee a converged 
solution* 
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without ILSS® 


Kith EI.PI 

A 

Simple 

Farturbatioa 

Finite Element 

Finite llesent 


BllipCc 

Solutioa 

First 

Final 

Final . . 


E@@ponoe 
(ref. 9) 


Iteration 

Result 

Result 


Simply Suppertcd ieasa Subjected to =• 

2.0 (Fq « 

1322 Ib/ira.) 

“ 1.0 

1.7852 

1.7854 

1.7852 

1.7856 

1.7682(3)^ 

* 2.0 

0.8472 

0.8660 

0.8621 

0.8460 

0.7108(4) 


1.6S57 

1.6503 

1.6563 

1.6S12 

1.5829(4) 

* 3.0 

1.4003 

1.4216 

1.4102 

1.3760 

1.2123(4) 


1.8217 

1.8314 

1.8226 

1.8002 

1.6743(4) 

dr 4.0 

1.8413 

1.8708 

1.84S3 

1.7846 

1.5871(6) 


2.1013 

2.1213 

2.0988 

2.0495 

1.8753(6) 

& 5.0 

2.260S 

2.2995 

2.252S 

2.1619 

1.9371(7) 


2.4361 

. 2.4673 

2.4236 

2.3432 

2.037(7) 


Clamped Subjected 

to P "1.0 
o 

(Fq " 3277 Ib/ln.) 

*1.0 

0.2118 

0.2165 

0.20% 

0.2091 

0.1772(3) 


1.4307 

1.4307 

1.4297 

1.4297 

1.4251(3) 

* 2.0 

0.8279 

0.8292 

0.8215 

0.8203 

0.7505(4) 


1.2987 

1.2990 

1.2942 

1.2936 

1.2743(4) 

± 3.0 

1 .0401 

1.0433 

1.0279 

1.0239 

0.9726(5) 


1.3232 

1.3248 

1.3127 

1.3099 

1.2694(3) 

* 4.0 

1.2183 

1.2247 

1.1979 

1.1S88 

1.1151(6) 


1.4101 

1.4142 

1.3910 

1.3836 

1.3197(6) 

* 5.0 

1.3938 

1.4042 

1.3619 

1.3457 

1.2513(8) 


1.5322 

1.5401 

1.5016 

1.4874 

1.4014(8) 


e. EffecCa of loisuitudinal defoEraeiosi aad ioercla. 

b. Number in brackees deuotes the auaber of iterations to get a converged 
oolutiou. 



closeness bet^en the earlier finite element results without ELDI, the 
simple elliptic response and the perturbation solution. The present finite 
elanent results Indicate clearly that the ELDI are to reduce the nonlinear- 
ity. 

Beams with various boundary conditions, including movable axial ends, 
are given in reference 8. Results obtained will be presented at the Second 
International Conference on Recent /yvances in Structural Dynamics, to be 
held April 9-13, 1983, at the Institute of iiound Vibration Research, South- 
hampton, England. 
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